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Abstract 

Boolean networks have been the object of much attention, especially since S. KaufF- 
man proposed them in the 1960's as models for gene regulatory networks. These systems 
are characterized by being defined on a Boolean state space and by simultaneous up- 
dating at discrete time steps. Of particular importance for biological applications are 
networks in which the indegree for each variable is bounded by a fixed constant, as was 
stressed by Kauffman in his original papers. 

An important question is which conditions on the network topology can rule out 
exponentially long periodic orbits in the system. In this paper, we consider systems with 
positive feedback interconnections among all variables (known as cooperative systems), 
which in a continuous setting guarantees a very stable dynamics. We show that for 
an arbitrary constant < c < 2 and sufficiently large n there exist n-dimensional 
cooperative Boolean networks in which both the indegree and outdegree of each variable 
is bounded by two, and which nevertheless contain periodic orbits of length at least c". 
In Part II of this paper we will prove an inverse result showing that any system with 
such a dynamic behavior must in a sense be similar to the example described. 

Keywords: Boolean networks, monotone systems, periodic solutions, mathematical biol- 
ogy, gene regulatory networks 
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The concept of a Boolean network was originally proposed in the late 1960's by Stuart 
Kauffman to model gene regulatory behavior at the cell level [12^ [T3] . This type of mod- 
eling can sometimes capture the general dynamics of continuous systems in a simplified 
framework, e.g. without the choice of specific nonlinearities or parameter values; see for 
instance pp. Boolean networks are known and used in several other disciplines such as 
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electrical engineering, computer science, and control theory, and analogous definitions are 
known under various names such as sequential dynamical systems [7] or Boolean difference 
equations [B]. 

An important class of continuous dynamical systems is that of so-called monotone sys- 
tems, which can be roughly characterized by the absence of negative feedback interactions 
[21 [15]. A special case is that of cooperative systems, in which there are no direct inhibitory 
interactions between any two variables. Monotone and cooperative systems have been used 
as a modeling tool for gene regulatory systems, for instance in [3]. The assumption of 
monotonicity is a stringent condition which ensures that the system behavior is remarkably 
stable: for instance, under mild additional assumptions the generic solution of a monotone 
dynamical system must converge towards an equilibrium. 

In the Boolean case, the class of cooperative systems can be described as that corre- 
sponding to maps that can be expressed using only AND and OR gates, i.e. with no use 
of negations. This can be easily seen by considering the disjunctive normal form of the 
Boolean maps. 

An important question in the study of cooperative Boolean networks is whether some 
of the stability properties of continuous cooperative systems have analogues in the Boolean 
case. For instance, does the assumption of cooperativity by itself limit the length of the 
longest cycle in an n-dimensional Boolean system? It was shown recently through simula- 
tions that random Boolean systems tend to have shorter periodic cycles if they are cooper- 
ative, or even if they are close to cooperative in the sense of having few negative feedback 
interactions; see [16] . and also [9l[l9]. Nevertheless, a straighforward use of Sperner's the- 
orem shows that a cooperative n dimensional Boolean system can have a cycle of length 
close to 2" for large n, see [8] and more recently [11^ 117]. 

One would like to know which additional assumptions rule out exponentially long pe- 
riodic orbits in cooperative Boolean systems. In [11] suitable adaptations of the notion of 
strong cooperativity [15] to Boolean systems were found that limit the length of periodic or- 
bits to 2v'"iog"(i+o(i)) 

or even to n, the dimension of the system. In the present manuscript 
we follow up on this question by considering a different class of cooperative Boolean systems 
in which both the indegree and the outdegree of the associated digraph is bounded. 

We need some definitions. An n-dimensional Boolean dynamical system or Boolean 
network is a pair (n.,g), where 11 = {0, 1}" and g : 11 — > 11. A state s{t) at time t will be 
denoted by s{t) = [si{t), . . . , Sri{t)], or simply s = [si, . . . , s„] if time-dependency is ignored. 
We will have 



The cooperative order on 11 is the partial order relation defined by s < r iff Sj < rj for 
alH G {1, . . . , n}. The system is cooperative if s{t) < r{t) implies s{t + 1) < r(i+ 1). 

We associate a directed graph D with vertex set {1, . . . ,n} with the system. A pair 
< i,j > is in the arc set of D iff there exist states s,r G 11 such that Si < ri and Sk = r^ 
for all A; 7^ i with the property that {g{si))j < {g{ri))j. We will say that the system is 
bi-quadratic if both the indegree and the outdegree of all vertices in D is at most two. 
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Already in his 1969 papers |12t ll3j. Kauffman focused his attention on Boolean networks 
where every variable can only be directly affected by a fixed number K of other variables. 
In the digraph associated to the network, this corresponds to limiting the indegree of every 
node to (at most) K. This corresponds to empirical findings about actual gene regulatory 
networks which show that most genes are directly regulated by a small number of proteins 
in a scale- free manner [H |18] . Other studies of biochemical networks show that only very 
few nodes are involved in the regulation of other chemicals. Thus large subnetworks of 
most biochemical networks of interest will also have the property that the outdegree of each 
node is bounded by a small integer. Bi-quadratic Boolean networks satisfy both of these 
restrictions with K = 2. Random Boolean networks with K = 2 have been extensively 
studied and tend to have dynamics in the ordered regime, which is characterized, among 
other properties, by the absence of exponentially long attractors (see |14j for a review). 
Thus it becomes a natural question whether one can prove, for cooperative bi-quadratic 
Boolean networks, a subexponential bound on the length of their periodic orbits, or at least 
a bound of the form c" for some constant c < 2. The following theorem shows that this is 
not the case. 

Theorem 1 Let c < 2 he arbitrary. Then for some sufficiently large n there exists an 
n- dimensional, bi-quadratic cooperative Boolean network which contains a periodic orbit of 
length at least c" . Moreover the digraph D associated with this network is strongly connected. 

The last sentence of Theorem [1] is of interest in connection with the results in [11] . 
There, we define a local version Dg of D for every state s as follows: A pair < i,j > is 
in the arc set of Ds iff there exist a state r S 11 such that either Sj < rj while s^ = r^ 
for all k ^ i, and we have {g{si))j < {g{ri))j, or < Sj while Sk = r^ for all k ^ i, and 
we have {g{ri))j < {g{si))j. It is shown that if X is a periodic orbit of an n-dimensional 
cooperative Boolean system such that Ds is strongly connected for every s € X, then 
\X\ < n (Theorem 25 of [11]). 

The proof of Theorem [1] uses a construction similar to a small Turing machine operating 
on a long circular tape. In part II of this paper we will show that if c is sufficiently close 
to 2, then all n-dimensional bi-quadratic cooperative Boolean networks with periodic orbits 
of length > c^ must contain a relatively small subsystem that can be considered a Turing 
machine operating on one or more tapes that retain the values of all other variables. 

The remainder of this note is organized as follows: In Section [Tj we introduce the 
main idea of the construction, but without requiring the system to be cooperative and 
bi-quadratic. In Section [2] we show how to modify the construction so that the network will 
also be cooperative, bi-quadratic and will have a strongly connected digraph. 

1 A Simple Counting Model 

In this subsection we consider a conceptual model of a (not necessarily bi-quadratic or 
cooperative) Boolean network with periodic orbits of length 2^, for arbitrary > 0. 
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We also discuss the problems that are involved in constructing such a network under the 
restrictions of Theorem [TJ Consider the states si, ... sat, and the system defined by 

Si{t) := Si+i{t - 1), i = l,...,N -1, 

s,^{t) := 7(si(t - l),mode{t - 1)). ^ ' 

One can think of 7 on a conceptual level as a Turing machine operating on variables 
numbered i = 1, . . . , N whose values are written on a circular tape. The variable mode can 
have one of two possible values for every t, namely mode = rotate, and mode = switch, and 
the function 7 is defined by 

7(x, rotate) = x, , > 

'y{x, switch) = 1 — x. 

Thus while mode{t) = rotate, iterating this machine will cyclically rotate the values of 
si, . . . ,sn- Whenever mode = switch, the machine also will rotate the variable values, but 
it will invert them at the site sn- 

Now let us define the value of the variable mode, in such a way that this machine 
behaves like a counter in base two. Let us require that at the times t = 0,N, 2N, 3N, . . . , 
mode{t) = switch. For all other times t, define 

moae[t) . | ^^^^^^^ if si(t-l) = 0. ^ ^ 

Thus the model turns into switch mode exactly at the times t = 0, N,2N, . . ., and it 
only returns back to rotate mode after Si(ti) = for some ti > t. The following lemma 
shows in what way this machine is a counter: if the states of the system encode numbers 
in binary format appropriately, then iterations are equivalent to the addition of one unit 
modulo 2^. 

Lemma 2 Given any state s of the model, define a{s) := si2^ + 522-*^ + . . . + S]\f2^~^ . Then 
a{s{N)) = a{s{0)) + 1 mod 2^. 

Proof: Consider an initial state s(0) and let j > be such that Sj(r/) = 1, for 1 < rj < j < 
N, and Sj+i(0) = 0. Note that a{s{0)) < 2^ — 1 in this case. We have mode{0) = switch 
by the definition above By ([2]), si(r/) = 1 for < r] < j — 1, si{j) = 0. Therefore 
mode{rj) = switch, for 1 < r] < j, and mode{j + 1) = . . . = mode{N — 1) = rotate. At time 
t = N, the variable values have completed a full rotation and returned to their starting 
points, except that = for 1 < r/ < j, Sj^i = 1, and Sj+2, ... sat are unchanged. Clearly 
a{s{N)) = a{s{0)) + 1 in this case. 

It remains to show the result for the case j = N, i.e. Sj(0) = 1, for every i = 1, . . . ,N. 
In that case mode{0) = mode{l) = . . . mode{N — 1) = switch by ([2]) and ([4]). In this way 
every value of the system is inverted at si from 1 to 0, so that Si{N) = for i = 1 . . . N. 
Therefore a{s{N)) = = a(0) + 1 mod 2^. □ 
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Corollary 3 The network given by equations (0), 0^, contains a periodic cycle of 
length at least 2^ . 

Proof: Since the variable mode is reset to switch for t = 0, N,2N, . . ., Lemma [2] apphes 
at each of these time points. Therefore one can start with s(0) = 0, and apply Lemma [2] 
successively to reach states s(0), s{N), s{2N), . . . , s{{2^ — 1)N), which are all different from 
each other. □ 

Importantly, the function 7 negates the values of the input x in switching mode. This 
appears to be an essential non-monotonic component (or negative feedback) of this system. 
Nevertheless, it is shown below that in fact one can rewrite our system in such a way that 
the resulting system is cooperative. 



1.1 A Generalized Counter 

Before proceeding with the proof of the main result, consider the following generalization 
of the simple counter above. Instead of individual Boolean values, each variable is now 
considered to be a vector with I > 1 Boolean entries, Sj = {s[, . . . , sj). We will treat Sj as a 
binary code for a nonnegative integer < 2'. At each time t, the system continues to be in 
one of two modes mode{t) = switch or mode[t) = rotate, but the function 7 is now replaced 
with a vector function T which we describe in the next paragraph. 

As before, when mode = rotate we let r(x, mode) := x. When mode = switch, and given 
X = {x\x^~^ . . . ,x^) 7^ (1, . . . , 1), let j be such that x^ = 1 for 1 < t] < j < I, x^~^^ = 0. 
Define y by letting := for 1 < 1] < j , letting y^'^^ := 1, and y^' := x^ for j + 1 < r] < I. 
Set T{x, switch) := y. If x = (1, . . . , 1), set T{x, switch) := (0, . . . , 0). In other words, the 
function T{x, switch) is defined as the addition of 1 to the vector x, in base 2 and modulo 2'. 

We define the generalized system 

Si{t) := Si+i{t - 1), i = l,...,N -1, 

S7v(t) := r(si(t - l),mode{t - 1)), ^ ' 

where F is defined as above. The variable mode{t) has the value switch for t = 0, A^, 2A^, . . . 
and for other values of t: 

mode{t) := | "^'ff " 7 1) = (1> • • • > 1)> (g) 

[ rotate, otherwise. 

Lemma 4 The network defined by equations contains a periodic cycle of length at 

least 2^K 

Proof: For {x\ ...,x^) G {0, 1}', define (5{x) := x^2^ + x'^2^ + . . . + x^^'^ . Note that 
[5 {T{x, switch)) = + 1 mod 2'. We follow an argument very analogous to Lemma [2] 
and Corollary El Let a{s) := /3(si)(2')0 + /3(s2)(2')^ + . . . + P{sn){2^)^-^ . Thus the vector 
. . . ,P{s]\f)) can be regarded as the representation of a{s) in base 2'. 
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As in the proof of Lemma [21 consider an initial state s(0), and let j > be such that 
Sri{0) = (1,...,1), for I < r] < j < N, and Sj+i(0) 7^ (1, . . . , 1). As before, we have 
mode{r]) = switch for < rj < j , and mode{j + 1) = ... = mode{N — 1) = rotate. At time 
t = N we have = (0, . . . , 0) for 1 < < j, as well as /3(sj+i) = /3(sj+i(0)) + 1, and 
Sj_|_2, . . . ,sj\f are unchanged from t = 0. Clearly a{s{N)) = a(s(0)) + 1. 

In the case that Sj(0) = (1, . . . , 1) for every i = 1, . . . ,N, it follows as before that 
mode(O) = mode(l) = ... = mode{N — 1) = switch. Therefore Si{N) = (0, ...,0) for 
i = l...iV, and a{s{N)) = 0. 

Repeating this process for s(0) = and t = N, 2N, . . . , as in Corollary [3l one finds 
states s of the system such that a{s) = 1,2,..., and which are therefore pairwise different. 
When Si = (1, . . . , 1) for all i, that is, when a{s{t)) = (2')^ — 1, this process reverts to 
a{s{t + N)) = 0. □ 



2 A Cooperative Counter 

In this section we carry out a construction which is analogous to that in Section [H but 
in which the underlying Boolean network is cooperative, bi-quadratic, and has a strongly 
connected digraph. We will need to define some auxiliary Boolean networks with designated 
input and output variables. 

Throughout this section let L > be an arbitrary even number, and consider the 
set A := {{ri, . . .ri) G {0, 1}^ | si + . . . + sl = L/2}. Define the special sequences 
START = (1, . . . , 1, 0, . . . , 0), i.e. L/2 ones followed by L/2 zeros, and similarly ACTIVE = 
(0,...,0, !,...,!). 

Lemma 5 Let g : A ^ A be an arbitrary function. There exists a Boolean network B with 
input vectors a = (ai, . . . ^ol), d = {di^d2), and output vector c = (ci, . . . ,cl), such that 
for some fixed m > the following equation holds for every t and a{t) € A, regardless of 
the initial condition of B : 

I 9{a{t)), ifd[t) = (1,0). 

Furthermore, the network B is cooperative, every node of its associated digraph has in- 
and outdegree of at most 2, and the indegree (outdegree) of every designated input (output) 
variable is zero. 

Proof: Define the set A := ^x{(0,l), (1,0)}, and the function G : A AhyGix, (1,0)) := 
g{x)^ G{x, (0, 1)) := x, for arbitrary x £ A. Since A is an unordered set, G can be extended 
to a cooperative function G : {0,1}^"*"^ {0,1}'^; see [H]. The result will follow from 
building a Boolean network that computes the function G. 

Consider a fixed component Gi : {0, 1}''""'"^ — > {0, 1} of G. By the cooperativity of 
this function, one can write it in the normal form Gj(yi, . . . ,yL+2) = ^i(yi, • • • ,yL+2) V 
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. . . V (yi, . . . , yL+2)) where each is the conjunction of a number of variables, i.e. 

. . . , yL+2) = Van A ... A Udji- This suggests a way of computing Gi: define Boolean 
variables ^j(t) := - 1)), and then let Gi{t) := V'Kt - 1) V . . . V ^{^{t - 1). Repeating 

this procedure for all components of G yields a Boolean network which computes G in 
m = 2 steps, and which is cooperative and has indegree (outdegree) zero for every input 
(output). 

In order to satisfy the condition that every node have in- and outdegree of at most two, 
we need to modify this construction by introducing additional variables. First, note that 
the outdegree of every input yi can be very large. One can define two additional variables 
which simply copy the value of yi{t), then four variables that copy the value of the previous 
two, etc. This procedure is repeated for each yi so that at least as many copies of each 
variable are present as appear in the expressions of all ipj. A similar cascade can be used to 
define each jpj and Gj so that each indegree is at most two. If ipj = yai f\ya2 ^Va^, say, then 
one can define zi{t) := yQ,^(t-l), Z2{t) := ya^it-l) Aya^it-l), ipl{t) := zi{t-l) Az2it-1). 
Similarly for longer disjunctions and each ipj and also similarly for Gi, in which case A is 
replaced by V at each step. This produces a computation of Gi in rrii steps for each i. 
Finally, after introducing further additional variables at each component i if necessary to 
compensate for unequal lengths of the expressions for ^pj, the Boolean vector G{yi, . . . , yL+2) 
can be computed in exactly m = max(mi, . . . ,7111) steps. □ 

Remark: Without loss of generality, we can assume that for every state variable s in the 
network B, there exists some input variable di or and a directed path from this input 
towards s. This is because if that wasn't the case, one could delete s from the system 
without altering equation ([7|). Cooperativity of G implies that G(0, . . . , 0) = (0, . . . , 0) and 
G{1, . . . , 1) = (1, . . . , 1) [11]; therefore each Gi is non-constant and no output variable will 
be deleted. Similarly, it will be assumed that for every state variable s, there exists an 
output variable Cj such that there is a directed path from s to q. 

Lemma [5] can be used to compute a function g which will be used in a way analogous to 
7 in equation ([2]). Similarly, we need to construct a 'switch' to determine when to turn the 
system into rotate mode, which is provided by Lemma [6] below. Note that Lemma [5] cannot 
be used for this purpose because the desired output depends not only on the current state 
of the input p{s) but on the whole history (of unknown length) of the input sequence since 
the last time when p{s) took the value START. 

Lemma 6 There exists fi > and a Boolean network D with input vector p = (pi, . . . ,pl), 
and output vector q = {qi,q2), such that the following holds for any initial condition of D. 
Consider any sequence of inputs p(0),p(l), . . . ,p{M), M > 1, such that 

i) p{s) e A, forO<s<M, 

ii) p{0) = START, and 

Hi) p{s) / START, for 0<s<M. 
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Let j >0 be such that p{s) = ACTIVE for l<s<j, p{j + 1) / ACTIVE (or p{l) 
p{M) = ACTIVE and j = M). Then 



Qis) 



(1,0), n<s<fj. + j, 
(0,1), i^i + j<s<fi + M, 



(8) 



Furthermore, the network B is cooperative, every node of its associated digraph has in- 
and outdegree of at most 2, and the indegree (outdegree) of every designated input (output) 
variable is zero. 



Pi 



U 



V 




Figure 1: The digraph of the network D which is used to compute the output qi from the 
input p. The formulas for each interaction (i.e. A, V) as well as the dependencies of U2 on 
u and V2 on v are omitted in this figure. 



Proof: The idea for this proof is the simple system c{t) = u{t — 1) V d{t — 1), d{t) = 
v{t — 1) A c(t — 1), with inputs u, v. This switch is turned on by letting both inputs u = 1 
and V = 1 for a short time, after which u can be turned to while v is left equal to 1. After 
letting V = for a short time, the switch resets and doesn't restart even if = 1 again. 

Let t = without loss of generality, the more general case being completely analogous. 
For the sake of clarity assume for now that < j < M, but the same construction allows 
for j = and j = M as described below. See Figure [1] which displays the circuit described 
below. Define for the moment u{t) := pi(t — 1) A . . . A PL/2{t — 1), v{t) := PL/2+ii't — 1) A 
. . .Apiit — l) (a modification of this definition with additional variables and indegree two is 
displayed in the figure and described below). Thus u{s) = 1 if and only if p(s — 1) = START, 
and v{s) = 1 if and only if p{s — 1) = ACTIVE, since by assumption p{s) G A. 

Define 

ui{t) := u{t - 1), U2{t) := u{t - 1) V ni(t - 1), Wi{t) := U2{t - 1), u^it) := u^{t - 1), 

Vl{t)=v{t-l), V2{t)=v{t-l)hVi{t-l), Vsit) ■=V2{t-l), Vi{t) ■.= Vs{t-l)VWl{t-l), 
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Wo{t) := ui{t - 1), wi{t) := Wo{t - 1) Avi{t - I), 



c{t) := Ui{t - 1) V d{t - 1), d{t) := v^it - 1) A c{t - 1). 

(Intuitively, U4 is a time-transposed copy of u where every 1 has been doubled due to 
the feed-forward loop at U2- Also, U4 is similar to a time-transposed copy of v where every 
has been doubled - the auxiliary variables Wi only play a role at a single time step as 
described below. The loop c ^ d forms the core of the switch in the system.) 

A simple calculation shows that tt4(4) = 7x4(6) = 1, U4{s) = for 5 < s < M -|- 4. 
On the other hand, since v{l) = 0,v{2) = . . . = v{l + j) = l,v{2 + j) = 0, we infer that 
V2{2) = 172(3) = 0, U2(s) = 1 for 3 < s < 2 j, ^2(3 + j) = '^2(4 + j) = 0. It follows that 
wi{3) = (since vi{2) = 0), and that i(7i(4) = 1 if and only if vi{3) = 1 (since wo{3) = 1). 
This in turn holds since j > 0. Also, wi{s) = for s > 4. 

We use the data for wi and ^3 to compute the values of V4. From wi{3) = ^3(3) = 0, it 
follows that ^4(4) = 0. From wi{4:) = 1 it follows that ^4(5) = 1, and using ^3 we similarly 
infer that ^4(5) = 1 for s = 4 < s < 4 -|- j. Also, ^4(5 -|- j) = 1)4(6 -|- j) = 0. 

We conclude that c(5) = 1, d(5) = 0, regardless of the values of c, d at earlier time steps. 
Since j > 0, one has c(6) = 1, d(6) = 1, and in general c{s) = d{s) = 1 for 5 < s < 5 -|- j. 
Then c(6 + j) = 1, d{6 + j) = 0, c{s) = d{s) = 0, ior 7 + j < s < 5 + M, and d(6 + M) = 0. 

In particular d{s) = 1 for exactly j time steps, 5 < s < 5 + j, and then d{s) = for 
6 + i<s<6-|- M. Since we want the variable qi to be equal to 1 during exactly j + 1 time 
steps, we define the additional variables 

W2{t) := Wo{t - 1), W3{t) := W2{t - 1), qi{t) := w^it - 1) V d{t - 1). 

Calculating that 1^3(5) = 1, ^3(5) = 0for5<s<5-|- M, we have qi{s) = 1, 6 < s < 6 + j, 
and qi{s) = 0,6 + j <s<7 + M. 

In order to define the variable q2, it suffices to make a construction dual to the previous 
one (recall that simply negating qi is not permitted). That is, define u{t) := PL/2+ii't — 1) V 
. . . V pL{t — 1), and v{t) := pi(t — 1) V . . . V PL/2{'t — 1), in such a way that u{s) = if and 
only if j>(s — 1) = START, and v{s) = if and only if p(s — 1) = ACTIVE. Define variables 
ui,vi etc. similarly as above, except that every A in the function definition is replaced by V 
and vice versa. Then it will necessarily follow that q2 = ^qi on the interval 6 < s < 6 -|- M. 
Using the value /i = 6, equation ^ is satisfied. 

The case j = is very similar as above. In this case iui(4) = (instead of for j > 0), 
774(4) = 774(5) = 0, and therefore d{s) = on all 6 < s < M -|- 6. Thus qi{6) = 1, and qi = 
for larger values of s. In the case j = M, one can compute 174(5) = 1 for 5 < s < M -|- 5. 
This allows the variables c{s),d{s) to remain equal to 1 up to and including s = M + 5. 
Therefore qi{s) = 1 up to and including s = 6 -|- M. 

Notice that this system is cooperative, and that all in- and outdegree requirements are 
satisfied except for the indegree of the variables u, v, u, v. These terms can now be replaced 
in a routine manner by a cascade of variables (see Figure [1]), in such a way that u{s) = 1 
if and only if p{s — r) = START, etc., for some r > 1. This will increase the delay fi but 
leave the computations and the other properties of this system unchanged. □ 
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We are ready for the construction of the cooperative counter described in the introduc- 
tion. This Boolean network is designed to rephcate the behavior of the system described 
by equations ([S]), while ensuring its cooperativity. In order to do so, we let / > be 
arbitrary and L > be an even positive integer, which is large enough that there exists 
an injective function x '■ {0, 1}' A, where A is defined as above. The cooperative net- 
work will contain L-dimensional vectors = {rj, . . . ,rf"), with values in A, which will be 
considered as proxy for states Sj = x~^{i"i) of the system (l5|), ([6]). 

We require that x(l, •••,!) = ACTIVE, and that START Im(x) (see the definitions 
of START and ACTIVE above). This is possible if L is large enough so that (^^2) > 2'- 
We also let x(0, ...,0) = (1, 0, . . . , 1, 0), and x(0, 0, . . . , 0, 1) = (0, 1, 0, 1, . . . , 0, 1). Having 
defined we define g : A ^ A as g{r) := xi^ix~^ {'''))) 1 for G Qi''') = ''' for 

other r € A. The function T is defined as in Section [TJ In particular, (^(START) = START. 

Using the function g defined above, we consider the cooperative networks B and D from 
Lemmas El and m Recall that B (D) has variables a,d (p) which are specifically designated 
as inputs, a variable c (g) specifically designated as output, and a 'processing delay' m (/x). 
The cooperative network, which will be denoted by S, is defined by B and D, together with 
the equations 

ri{t) := ri+i(t - 1), i = m + 2,m + 3, . . . , N, , . 

rN+i{t):=c{t-l), 

and 

a{t) := rm+2{t - 1), 

d{t):=q{t-l), (10) 
p{t) := rm+i,+2{t)- 




Figure 2: The network interconnections of the system S given by B, D, and equations Q, 
(|10|) . The variables ri, . . . , r^+i are displayed in a box to indicate that they are not part 
of S but only included in the proof of Theorem [HI 

See Figure[2]for an illustration. Since both of the subnetworks used in the construction of 
this system contain only the Boolean operators A, V in their expression (and no negations), 
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it follows from ([9]) and (jlOp that the same is the case for the full network, hence the system 
is cooperative. 



Proposition 7 The digraph of the Boolean network S is strongly connected and bi-quadratic. 

Proof: The fact that every in- and outdegree is at most 2 follows directly from equations 
([9]) , ([To]) and Lemmas [5] and [H taking into account that the indegree (outdegree) of every 
input (output) variable is zero within their respective subnetwork. See also Figure [2j 

In order to show the strong connectivity of the digraph, first we show that there exists 
a directed path from every node in the network to the node qi, the first component in the 
output of -D. It is clear from the circuit defining D that every input variable pi has a path 
connecting to qi (the first L/2 components of p through the variables u,ui, . . . and the last 
L/2 components through v,vi,...). Therefore every component of every variable rj can 
reach qi as well. By the remark after Lemma O the same applies to every variable of c, and 
thus to every variable in the subnetwork B. Thus the same applies also to q2, and hence to 
every state in the subnetwork D. 

Now we show that there exists a path from qi to every node in the network. Suppose 
first that there exists cj such that neither di or d2 contains a path towards Cj. This would 
imply that gj{x) = xj for every argument x £ A, hy equation d?]). But we have 



which is a contradiction. Thus for every j, there exists a path from either di or d2 to Cj 
(and therefore from qi or q2 to Cj). 

Since there exists a path from qi to q2, it follows that there is a path from qi to every Cj. 
Thus every component of every state Vi, p, and a can be reached by a path from qi. Every 
state in B can be reached from di and hence qi, once again by the remark after Lemma [Sj 
the same applies to q2, and every state in the subnetwork D. □ 

Theorem 8 Let L > be an even number such that {j^i^ > 2'. Then system S has a 




periodic orbit of length greater than or equal to 2 . 

Proof: For the purposes of this proof, we extend the system with the auxiliary variables 
ri, . . . , Tm+i, defined by rj(t) := rj-|_i(t — 1), z = l,...,m + l; see Figure [2j These variables 
cannot change the length of the original system's periodic orbits (since they don't feed back 
into it), but they can nevertheless be used for the study of the network. 

Suppose that the system is initiated at time to, and let t > to + ^ + 2. Then ri{t — 1) = 
rm+2{t - m - 2) by dS]). But then a{t - m - 1) = rm+2{t - m - 2) = ri{t - 1), by (fin]) . 
By Lemma O c(t — 1) is equal to either ri{t — 1) or g{ri{t — 1)), depending on whether 
d{t — m — 1) = (0, 1) or (1,0) respectively. Since r7v+i(i) = c(t — 1), we have 



ff(i, 0, . . . , 1, 0) = x(r(o, . . . , 0)) = x(o, . . . , 0, 1) = (0, 1, 0, 1, . . . , 0, i) 





ri{t — 1), if mode{t — 1) = rotate, 
g{ri{t — 1)), if mode{t — 1) = switch. 



(11) 
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where the auxihary Boolean variable mode{t) is defined as mode{t) := switch if d{t — m) = 
(0, 1) and mode{t) := rotate if d{t—m) = (1, 0). The variable mode, similarly as ri, . . . , r^+i, 
is defined merely for the purposes of this proof, and it does not form part of the network 
itself. 

Suppose now that tg ^ — fi — 2. At time 0, assume that r^r+i = START, and 
/ START, for I < rj < N. Let j > be such that = ACTIVE for 1 < 77 < j < iV, 
and r^+i / ACTIVE. We show that 

mode(r]) = switch, < rj < j; mode{rj) = rotate, j + 1 < f] < N. (12) 

To see this, note that by Q rm+^+2(??) 7^ START, for —m — n—l<r]<N — m — fi — 2. 
Since (START) = START, it also follows that START = ri(0) = c(-l) = a{-m - 1) = 
i"m+2{—m — 2), and rm+/x+2(— — — 2) = START. Thus setting t = —m — fi — 1, one has 
p{t) = START, p(r?) / START for t < ry < t + iV, ^(77) = ACTIVE for t + 1 < r? < i + j, 
and j5(t + j + 1) / ACTIVE. Applying Lemma[6]with M := N, we have that q{ri) = (1, 0), 
for t + ^ < < t + /J, + j, and g(r?) = (0, 1), i + ;U + j < < t + /j, + iV. Equation ([T2D then 
follows directly from the definition of d and the mode variable. It is analogous to verify 
that ^ also holds in the case j = N, i.e. when r^(0) = ACTIVE for 1 < r? < iV. 

Note that using equation (112p we can fully calculate r{N + 1), namely r^{N + 1) = 
g{r^{{))), for 1 < 7? < j + 1, and r^(A^ + 1) = r^(0) for j + 1 < ?? < A^; also, necessarily 
rN+i{N + 1) = START regardless of mode{N), since ^(START) = START. The same 
process can be repeated starting at time A^ + 1, 2(A^ + 1), etc., since necessarily = START 
can still only hold for 77 = + 1. 

An appropriate initial condition to reach the above situation can be given as follows. 
Let tQ = -{N + 1), and let r^(to) = x"HO, • • • , 0), for 1 < < N. Let rAr+i(to) = START. 
Finally, let B {D) be initialized with m (/x) successive inputs of a = x~^(0) • • • ) 0)) = (0) 1) 
{p = X~"'^(0, • • • , 0)). This way for t = we guarantee that tat+i = START, r^(to) = 
X~"'^(0, . . . , 0) for 1 < ?7 < A^, and importantly, to ^ — fi — 2. 

Finally, under our standing hypotheses to ^ ~ A* ~ 2, rjv+i(0) = START, and 
r^(0) 7^ START for 1 < rj < N . Define the following initial conditions for the system ([5]), ([6]): 
Sri{0) := X'^ifriiO)), i = 1...N. After calculating j as before, < j < A, we have 
seen that s»y(A) = r(s,,(0)) for 1 < rj < j + 1, and Sr^{N) = s,,(0) otherwise. From the 
discussion above, it follows that x~^(?'r?(A + 1)) = s^(A) for 1 < r] < N . This equivalence 
between the two systems implies in particular that the states r{t) are pairwise different for 
t = 0,A + 1,2(A + 1),...,(2^^ -1)(A + 1). The result follows. □ 

2.1 Proof of Theorem [1] 

We can use Proposition [7] and Theorem [8] to prove the theorem stated in the introduction. 
Let < c < 2 be arbitrary. We prove first that there exist L > even and and integer I > 
such that 
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(l/2) > 2' > (1^) 

The second inequality is equivalent to L/l < In 2/ In c; thus let L = wl, for some fixed 
l<t(;<ln2/lnc (for large enough /, L can then be rounded up to the nearest even number 
while satisfying this inequality). Using Stirling's formula, we have (^/2) > l^phxL for 
large enough L, where < u < 1 is arbitrary and fixed. The first inequality in (|13|) is satisfied 
if V 2^ /^2ttL > 2'. But after replacing L = wl this is equivalent to 2^"'"-'^)' > v~^^/2-kwI. 
Clearly this inequality is satisfied for sufficiently large /, hence (jlSp follows. 

The first inequality is now used to carry out the construction of system S, which by 
Theorems [7] and [8] is cooperative and bi-quadratic with strongly connected digraph, and has 
a periodic orbit of length greater than or equal to 2^'. 

It remains to show that 2^' > c" for large > 0, where n is the dimension of the 
system. Let T be the total number of variables in the subnetworks D,B. Note that T 
depends only on L, I, and not on A^. Then n = (N + 1 — (m + 1))L + T = NL — rriL + T. 
Notice that c" < 2^' if and only if {NL - ml + T)\nc<Nl In 2, which holds if and only if 

mL — T 
L In c < / In 2 H — In c. 

But this equation is satisfied for large enough A^, since Line < nn2 by (jl3p . □ 



References 

[1] Albert R. and Othmer, H.G. (2003) The topology of the regulatory interactions predicts 
the expression pattern of the segment polarity genes in Drosophila melanogaster, J. 
Theor. Biol. 223:1-18. 

[2] Angeli, D. and Sontag, E.D. (2003) Monotone control systems, IEEE Trans. Autom. 
Control 48:1684-1698. 

[3] Angeli, D., Ferrell, J.E., and Sontag, E.D. (2004) Detection of multistability, bifurca- 
tions, and hysteresis in a large class of biological positive- feedback systems, Proc. Nat. 
Acad. Sci. 101:1822-1827. 

[4] Arnone, M. I. and Davidson, E. H. (1997) The hardwiring of development: organization 
and function of genomic regulatory systems, Develoment 124:1851-1864. 

[5] Colon-Reyes, O., Laubenbacher, R., and Pareigis, B. (2004) Boolean monomial dynam- 
ical systems, Ann. Comb., 8(4):425-439. 

[6] Dee, D. and Ghil, M. (1984) Boolean difference equations, I: formulation and dynamic 
behavior, SIAM J. Appl. Math. 44(1):111-126. 



13 



[7] Garcia, L., Jarrah, A.S., and Laubenbacher, R. (2006) Sequential dynamical systems 
over words, Appl. Math. Comput. 174:500-510. 



[8] Gilbert, E.N. (1954) Lattice theoretic properties of frontal switching functions, J. Math. 
Phys. 33:57-67. 

[9] Greil, F., and Drossel, B. (2007) Kauffman networks with threshold functions, Europ. 
Phys. J. B 57:109-113. 

[10] Hernandez Toledo, R.A. (2005) Linear finite dynamical systems. Communications in 
Algebra 33:2977-2989. 

[11] Just, W. and Enciso, G. (2007) Analogues of the Smale and Hirsch theorems for coop- 
erative Boolean and other discrete systems, preprint, arXiv:math/0711.0138vl. 

[12] Kauffman, S.A. (1969) Homeostasis and differentiation in random genetic control net- 
works. Nature 224:177-178. 

[13] Kauffman, S.A. (1969) Metabolic stability and epigenesis in randomly constructed 
genetic nets, J. Theor. Biol. 22(3):437-467. 

[14] Kauffman, S. A. (1993) Origins of Order: Self- Organization and Selection in Evolution. 
Oxford U Press, 1993. 

[15] Smith, H.L. (1995) Monotone dynamical systems, Math Surveys and Monographs, 
AMS, Providence, RI. 

[16] Sontag, E.D., Veliz-Cuba, A., Laubenbacher, R., and Jarrah, A.S. (2007) The 
effect of negative feedback loops on the dynamics of Boolean networks, preprint, 
[aFXiv:0707.346 8v2 [q-bio.QM] 

[17] Sontag, E.D. (2007) Monotone and near-monotone biochemical networks, Journal of 
Systems and Synthetic Biology, to appear. 

[18] A.H.Y. Tong et al. (2004) Global mapping of the yeast interaction network. Science 
303:808-813. 

[19] Tosic, P.T., and Agha, G. (2004) Characterizing configuration spaces of simple thresh- 
old cellular automata. Springer- Verlag LNCS series, vol. 3305, pp. 861 - 870. 



14 



